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1. INTRODUCTION 


This is a versatile VUORTRAN subroutine for performing two 
dimensional plotting operations on 4a line printer, The subroutine 
explicitly provides for X-Y plotting, histograms, frequency 
distributions, and cumulative frequency distributions. Plotting of 
specified functions, preparation of contingency table matrices from 
numerical codes, and several other operations are also available 
depending upon the user's knowledge of programming. Although several 
programs are available for each of these operations, we feel that 
there is an advantage in having a single subroutine which can perform 
all of these with only minor modifications in the control cards and 
calling procedures. This uniformity is especially advantageous for 
users with a limited knowledge of programming. 


All testing and example runs were performed under the 
Version II disk monitor system on an IBM 1130 computer with 8K, 
16-bit words of core storage, IBM 1442 card read punch, IBM 1132 
printer, and a single magnetic disk; however, the subroutine does not 
require the auxiliary storage capacity of the disk. 


The senior author is responsible for the systems analysis and 
initial programming, while the junior author did extensive debugging 
and testing of the program, 


Il. PURPOSE 


This subroutine subprogram is designed to make a line 
printer, attached as an output device to a computer, serve as a 
versatile two-dimensional plotter. This is possible because of the 
approximate correspondence between one vertical space and two 
horizontal spaces on such a printer. The horizontal scale is 
actually a little larger than the vertical scale, and this should be 
kept in mind when interpreting plots produced by this program. 


The program provides for the construction of X-Y plots, 
histograms, frequency distributions, cumulative frequency dis— 
tributions, and contingency table matrices for numerical data. It 
can also be used for plotting functions, The only restriction is 
that the horizontal scale must increase from left to right, and the 
vertical scale must increase from bottom to top. An observed value 
is plotted at the nearest intersection of grid lines. A limited 
control over the treatment of boundary values is provided via a 
control card parameter. The user is free to select the symbols to 
be plotted. For X-Y plotting, the program can easily be made to plot 
the number of observations falling at each intersection of grid lines 
(up to and including a frequency of 99). 


The horizontal axis of the graph is the first row of plotted 
characters above the column marks, and the vertical axis of the 
graph is the first double column of characters to the right of the 
line marks. Every fifth value on each scale is printed. Since 


these scale values are printed under 4 fixed point furmat, some 
scaling of the data may be necessary or desirable when dealing with 
very large or very small values. The grid lines are omitted from the 
graph to improve readability. 





The number of observations plotted on the graph is printed 
in the lower left-hand corner of the page containing the graph. The 
number of observations which could not be plotted because they were 
beyond the scale limits is printed (if any) on the page immediately 
following the graph. Two sizes of graphs are available. The 
smaller size, having a resolution of | in 40, fits on an 84 X 11-inch 
page when the margins are trimmed. The larger size, having a 
resolution of 1 in 50, is 10 by 12 inches in size. An option is 
available for multiple copies. 


III. PROCEDURE 
The user must supply the following items: 


A. A short FORTRAN main program which specifies the form 
of the input data and invokes the necessary calls to 
SUBROUTINE PLOT. 


B. The monitor cards necessary for compilation and 
execution of the main program and subprogram. 


C. Seven control cards. The first three control cards 
are used to specify the desired plotting symbols. 
For convenient reference, we refer to these as the 
SYMBOLA, SYMBOLB, and SYMBOLC control cards, respec- 
tively. The fourth, or XTITLE, control card contains 
title information to be printed below the horizontal 
axis of the graph. The fifth, or YTITLE, control card 
contains title information to be printed at the top of 
the graph. The sixth, or OPTION, control card is used 
to set program parameters. These first six control 
cards are placed in front of the data deck. The seventh 
control card, called the STOP card, is placed immediately 
behind the data deck. It contains a flag number which 
indicates the end of the input data for the current run. 


D, The deck of data cards. 
An assembled input deck is shown schematically in Figure 1. 


The main program, monitor cards, and control cards are explained 
more fully in subsequent sections. 


Moni ter 


Cords 





‘4/4 TOB 


Figure 1. Diagram of input deck, including monitor cards for the 
IBM 1130 disk monitor system. Broken outline indicates 
that item should be omitted if SUBROUTINE PLOT has been 
stored previously on magnetic disk. 


ty 
IV. MAIN PROGRAM 


The general form of the mainline program for X-Y plots, 
histograms, frequency distributions, cumulative frequency distri- 
butions, and contingency table matrices is shown in Figure 2. In 
all cases, the user must alter the FORMAT statement to suit his 
input data. Also note the number 2 in column 12 of the READ state- 
ment. This is appropriate for an IBM 1130 computer having an IBM 
1442 card read punch. If the computing system is not of this type, 
replace the number 2 with the appropriate reference number for the 
card reader. 


For X-Y plotting amd constructing the matrices of contin- 
gency tables, the letter X in column 17 of the READ statement refers 
to the variable to be plotted on the horizontal axis, and the letter 
Y in column 22 refers to the variable to be plotted on the vertical 
axis. If the X-value for a given observation precedes the Y-value 
on the data card(s), it is only necessary to revise the F@RMAT card. 
If the Y-value precedes the X-value, it is also necessary to inter- 
change the positions of X and Y in the READ statement. 


For histograms, frequency distributions, and cumulative 
frequency distributions, the "Y(I)," portion of the READ statement 
is irrelevant. If the READ statement is left as shown, a dummy 
field description must be provided for Y in the FORMAT statement. 
Alternatively, the user can delete "Y(I)," from the READ statement 
and the corresponding field description from the FORMAT statement. 


The mainline program for plotting functions is somewhat 
different. The general form for this purpose is shown in Figure 3. 
The user must rewrite the statement which begins Y(I)= so that it 
defines Y(I) in terms of X(I) for the particular function being 
plotted. The corresponding statement in Figure 3 is included only 
for illustrative purposes and is the version that is appropriate to 
example #5, Obviously, this use depends on the ability of the user 
to express the desired function as a FORTRAN program statement. 


When used with an IBM 1130 computer having 8K words of 
core storage, the main program must necessarily be short since 
SUBR#UTINE PLOT requires most of the available storage capacity. 
When more core storage is available, the main programs shown can be 
incorporated into larger programs. The CALL LINK feature of the 
IBM 1130 disk monitor system can also be used to extend the size 
of the main program; however, the COMMON area available for passing 
information will be severely limited for 8K systems. 


Col. 7 


DIMENSION A(18) ,B(18) ,X(20) ,Y(20) ,IA(50,50) ,1B(101) ,1C(10) 
Ic(1)=1 
GO Th 4 
1 IF(IC(7)-1)2,5,2 
2 NI=IC(2) 
READ (2,3) (X(I) ,Y(I) ,I=1,NI) 
3 FORMAT(F1O.4,10X,F10.4) 
4 CALL PLOT(A,B,X,Y,C,D,E,F,G,H, IA, 1B, IC) 
Ic(1)=2 
CH TH 1 
5 CALL EXIT 
END 


Figure 2. Generalized main program for X-Y plots, histograms, frequency 
distributions, and contingency table matrices. Monitor cards 
are not included. 


Col. 7 
‘ 


DIMENSION A(18) ,B(18) ,X(20) ,¥(20) ,IA(50,50) ,IB(101) ,1¢(10) 
Ic(1)=1 
N=0 
1 CALL PLOT(A,B,X,Y,C,D,E,F,G,H,IA,IB, IC) 
Ic(1)=2 
IF(IC(7)-1)2,6,2 
2 Dd 3 I=1,10 
X(L)=CHI*D 
Y(I)=.005454154*X (1) *X(I) 
3 N=N+L 
IF(N-1C(5))5,5,4 
x(1)=¢ 
GO TO 1 
CALL EXIT ' 
END 


Dur 


Figure 3. General form of main program for plotting functions. 
Monitor cards are not included. 
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V. CONTROL CARDS 
The general purposes and positioning of the control cards 


are described in section III and Figure 1. Detailed instructions 
for completing the control cards are given in Table 1. 


Table 1. 
Card name Columns Information 
SYMBOLA 1 through 2 This set of two symbols is 


used to plot distribution curves 
and functions and to mark the 
upper edge of the bars in a 
histogram. For simple X-Y plot- 
ting and contingency table 
matrices, this set of symbols is 
plotted at any grid intersection 
containing a single observation. 


cy 3 through 4 This set of two symbols is 
used to shade areas of histogram 
bars and the area beneath 
distribution curves. If these 
two columns are blank, there will 
be no shading. For simple X-Y 
plotting and contingency table 
matrices, this set of symbols is 
printed at any intersection con- 
taining two observations. For 
plotting functions, these two 
columns should be blank. 


my 5 through 80 For histograms, distribu- 
tions, and functions, these 
columns should be blank. For 
simple X-Y plotting and contin- 
gency table matrices, the set of 
symbols in columns 5 and 6 is 
plotted at any grid intersection 
having three observations, the 
set of symbols in columns 7 and 
8 is plotted at any intersection 
having four observations, ..., the 
set of symbols in columns 79 and 
80 is plotted at any intersection 
having 40 observations. 


1 through 80 


1 through 38 


39 through 40 


-8- 
(continuation) 
Information 


For histograms, distri- 
butions, and functions, the 
entire SYMBOLB card should be 
blank. However, the blank 
card cannot be omitted. For 
simple X-Y plotting and con- 
tingency table matrices, the 
set of symbols in colums 1 
and 2 is plotted at any grid 
intersection having 41 observa— 
tions, the set of symbols in 
columns 3 and 4 is plotted at 
any intersection having 42 
observations, ..., the set of 
symbols in columns 79 and 80 
is plotted at any intersection 
having 80 observations. 


For histograms, distri- 
butions, and functions, the 
entire SYMBOLC card should be 
blank. Again, the blank card 
cannot be omitted. For simple 
X-Y plotting and contingency 
table matrices, the set of 
symbols in columns 1 and 2 is 
plotted at any grid inter- 
section having 81 observations, 
the set of symbols in columns 
3 and 4 is plotted at any 
intersection having 82 observa- 
tions, ..., the set of symbols 
in columns 37 and 38 is plotted 
at any intersection having 99 
observations. 


These two columns should 
always be blank because this 
is the set of symbols plotted 
at any grid intersection con- 
taining no observations. 


fable 1. (continuation) 
Card name Columns Information 
SYMBOLC 41 through 42 Yor histograms, distri- 


butions, and functions, these 
two columns should be blank, 
For simple X-Y plotting and 
contingency table matrices, 
this set of symbols is plotted 
at any intersection having 
more than 99 observations, 


si 43 through 80 Any information in these 
columns is ignored by the com- 

puter. 
XTITLE 1 through 72 The information in these 


columns is reproduced starting in 
column 21 at the bottom of the 
graph. Any information on the 
remainder of the card is ignored. 


YTITLE 1 through 72 The information in these 
columns is reproduced in the 
corresponding columns at the top 
of the graph. Any information 
on the remainder of the card is 
ignored, 


OPTION 1 through 10 For all purposes except the 
plotting of cumulative distri- 
butions, this is the midpoint of 
the lowest scale interval on the 
horizontal axis. For cumulative 
distributions, this is the upper 
limit of the lowest scale interval. 
This number should be right justi- 
fied, and the decimal point 
should be punched. If, however, 
the decimal point is omitted, 
columns 7 through 10 will be 
interpreted as decimal places. 


m il through 20 The size (width) of the scale 
interval on the horizontal axis 
should be specified in these col- 
umns. This number should be right 
justified, and the decimal point 
should be punched. If, however, 
the decimal is omitted, columns 
17 through 20 will be interpreted 
as decimal places. 


Card name 


OPTION 


Table L. 
Columns 


21 through 30 


31 through 40 


41 through 50 


53 


-10 
(continuation) 


In all cases, this is the 
midpoint of the lowest scale 
interval on the vertical axis. 
This number should be right 
justified, and the decimal 
point should be punched. in 
the absence of a decimal point, 
columns 27 through 30 will be 
interpreted as decimal places. 


The size of the scale inter~ 
val on the vertical axis should 
be specified in these columns. 
This number should be right justi- 
fied, and the decimal point should 
be punched. In the absence of a 
decimal point, columns 37 through 
40 will be interpreted as decimal 
places. 


Place here a number which is 
not found among the values to be 
plotted on the horizontal axis. 
This number should be right 
justified and the decimal point 
should be punched. In the 
absence of a decimal point, 
columns 47 through 50 will be 
interpreted as decimal places. 
This number must appear again on 
the STOP control card (at the end 
of the data deck). It serves as 
a flag to indicate when the last 
data card has been read. The only 
other restriction on its value is 
that the number of digits plus the 
decimal point cannot exceed the 
number of columns allotted to a 
field of the variable to be 
plotted on the horizontal axis. 


Normally this column is left 
blank, in which case any given 
plotting interval will include 
its upper boundary but not its 
lower boundary. If the number 1 


ati 


Table 1. (continuation) 


Card name Columns Information 
OPTION 53 (cont'd) is placed in this column, the 


lower boundary will be in- 
cluded in the interval but not 
the upper boundary. 


a 55 through 56 This is an integer num- 
ber (right justified with no 
decimal) specifying the num- 
ber of observations per data 
card. For simple X-Y plot- 
ting and contingency table 
matrices, a set of horizontal 
and vertical coordinates is 
counted as a single observa- 
tion. This number must be 
dess than or equal to 20 . 
For plotting functions, this 
column can be blank, 


» 59 If the purpose is to 
construct a histogram or 
simple frequency distribution, 
piace the number | in this 
column. If the purpose is to 
construct a cumulative fre- 
quency distribution, place the 
number 2 in this column. For 
any other purpose, place the 
number 3 in this colum. 


" 62 This is the device code 
number (symbolic unit number 
or data set reference number) 
assigned to the printer on 
which the graph is to appear. 
The code number for the IBM 1132 
printer is3 . 


65 This column is used to 
specify the size of the graph. 
If the number 1 is placed in this 
column, the graph will have a 
resolution of 1 in 40 and will 
fit on an 84 X Ll-inch page when 
the margins are trimmed. If the 
number 2 is placed in this 
column, the graph will have a 


=i 3. 


Table 1. (conclusion) 


Card name Columns Intormation 
OPTION 65 (cont'd) resolution of 1 in 50 and 


will be approximately 10 by 12 
inches in size. 


vs 68 This column is used to 
specify the number of copies 
desired, including the orig- 
inal. If this column is left 
blank, it will be assumed that 
a single copy is desired. 


STOP varies The number placed in 
columns 41 through 50 of the 
OPTION card is repeated on this 
last control card, but not 
necessarily in the same columns, 
Its positioning here depends on 
the positioning of information 
on the data cards. The number 
should be positioned on this 
card exactly as if it was the 
horizontal coordinate of an 
extra observation, With 
several observations per data 
card, the STOP card and the 
last data card will often co~ 
incide. As mentioned above, 
this number serves as a flag to 
indicate the end of the data 
deck. 


Any questions regarding preparation or use of control cards 
should be clarified by reference to the examples. 


VI. IBM 1130 DISK MONITOR CARDS 


The monitor cards necessary for processing this program under 
the IBM 1130 disk monitor system are indicated in Figure | and Figure 4. 
The positions of these cards in the input deck are shown in Figure 1, 
and their formats are shown in Figure 4. For a different monitor 
system, the appropriate manuals should be consulted. 


Column 
number 


NTCWONTZTOWMONSTCOBDONSTOWDONTOCBDONTONMDON 
SCOSCCH AAA ANNNNNMMNMNTSTSTST TST NNN HNHH CO 
// I@B 


*EXTENDED PRECISION 

*@NE WORD INTEGERS 

// DUP 

*STORE WS UA PLOT 
*IOCS(CARD,1132 PRINTER) 
*LIST SOURCE PROGRAM 


// XEQ 


Figure 4, Positioning of information on IBM 1130 disk monitor cards. 
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VIL. SUBRPUTINE PLAT SOURCE PROGRAM 


The meanings of the FORTRAN labels used in SUBR@UTINE PLOT 
are explained in Table 2. 


Table 2. 
FORTRAN label Meaning 
DEL This is the difference between 


XMIN and the upper limit of the 
first scale interval on the hori- 
zontal axis. It determines whether 
the plotting will be over midpoints 
or upper limits of scale intervals 
on this axis. 


FL This is working storage. 
I This is an iteration variable. 
Ic This array is used to pass 


program parameters between main 
program and subroutine. 


ICHAR This array contains the 101 
sets of two characters used in 
plotting. 

IDIF This is working storage. 

IDROP The value 0 means include 


upper boundary but not lower in 
scale interval. The value 1 means 
include lower boundary but not upper 
in scale interval. 


1G The value 1 means enter the sub- 
routine at statement number 5000. The 
value 2 means enter the subroutine at 
statement number 50 . 


IGRID This is the 50 x 50 array rep— 
resenting the intersections of grid 
lines. 


Ip This is the reference number of 
the printer on which output is to 
appear. 


IQ This is working storage. 


FORTRAN Label 


1s 


ISIZE 


ISTOP 


ITP 


ITYPE 


1Z 


KQUNT 


Lgcx 


Table 2. 


15. 
(continuation) 


Meaning 


This is working storage 


The value | for this variable 
means smaller version of graph. 
The value 2 for this variable means 
larger version of graph. Subsequent- 
ly, ISIZE indicates the number of 
scale intervals on an axis. 


This variable is set equal to 1 
when the reading of data is complete. 


This variable is used only for 
frequency distributions and histo- 
grams. It indicates the highest 
frequency which can be represented 
on the vertical axis. 


The value 1 for this variable 
means operation is construction of 
histogram or simple frequency dis- 
tribution. The value 2 for this 
variable means operation is con- 
struction of a cumulative frequency 
distribution. The value 3 for this 
variable means operation is anything 
not covered by the previous two 
codes. 


This is the subscript index of 
the zero point on the vertical scale 
of histograms and frequency distribu- 
tions. 

This is an iteration variable. 


This is an iteration variable. 


This is the number of observa- 
tions plotted on the graph. 


This is an iteration variable. 
This is the horizontal coordin- 


ate of an observation in terms of 
scale intervals. 


FORTRAN label 


LcCY 


NCOPY 


NI 


SAVE 


NgUT 


NPRNT 
SAV 


STOB 


XMIN 


XSTEP 


XTITL 


Table 


2 


=L6= 
(continuation) 


Meaning 


This is the vertical coordin- 
ate of an observation in terms of 
scale intervals. 


This is the reference number 
of the card reader. 


This is the number of copies 
of the graph to be printed. 


This is the number of observa- 
tions per data card. 


This is working storage. 


This is the number of observa- 
tions falling outside the scale 
limits. 


This is an iteration variable. 
This is working storage. 


This is the flag number which 
indicates the end of the data 
deck. 


This is a buffer array which 
holds the horizontal coordinates of 
NI points while they are being 
processed, 


This is either the midpoint or 
the upper limit of the lowest scale 
interval on the horizontal axis, 
depending upon the purpose. 


This is the scale interval on 
the horizontal axis. 


This array contains title 
information to be printed at the 
bottom of the graph. 


This is a buffer array which 
holds the vertical coordinates of 
NI points while they are being 
processed. 


=17= 


Table 2. (conclusion) 


FORTRAN Label Meaning 

YMIN This is the midpoint of the 
lowest scale interval on the verti- 
cal axis. 

YSTEP This is the scale interval on 
the vertical axis. 

YTITL This is an array containing 
title information for the top of 
the graph. 


The FORTRAN source program for SUBR@UTINE PLOT appears on 
the pages immediately following. The position of column 7 is marked 
at the top of each page. The cards are numbered sequentially in 
columns 78 through 80. 


Col. 7 
' 


SUBROUTINE PLOT(YTITL,XTITL,X,Y,XMIN,XSTEP ,YMIN,YSTEP ,STOB,DEL, 
1IGRID, ICHAR, IC) 
DIMENSIQN IGRID(50,50) ,ICHAR(1O1) ,YTITL(18) ,XTITL(18) ,X(20) ,¥(20), 
11¢(10) 
IGG=1C(1) 
NI=IC(2) 
ITYPE=IC(3) 
IG=1C(4) 
ISIZE=IC(5) 
NCOPY=IC(6) 
ISTOP=1C(7) 
IDROP=1C (8) 
NUT=1C (9) 
K@UNT=IC (10) 
C ....DETERMINE ENTRY POINT. 
GP T4(5000,50) ,1Go 
5000 M=2 
READ (M,1) ICHAR, XTITL, YTITL,XMIN,XSTEP , YMIN,YSTEP ,ST@B, IDROP,NI, 
LITYPE,19,ISIZE ,NCOPY 
1 FORMAT (40A2/40A2/21A2/,2(18A4,/) ,5F10.4,6(1X,12)) 
C ....ADJUST GRID T@ SPECIFIED SIZE. 
GP TO (601,602) ,ISIZE 
601 ISIZE=40 
GO TO 603 
602 ISIZE=50 
C ...,ASSUME SINGLE COPY @F GRAPH IF N@T SPECIFIED OTHERWISE. 
603 IF (NCOPY)610,609,610 
609 NCOPY=1 
C ... EXCEPT FOR CUMULATIVE, PL@T OVER MIDP@INTS @N THE ABSCISSA. 
610 IF(ITYPE-2)102,101,102 


ool 
002 
003 
004 
005 
006 
007 
008 
009 
010 
oll 
012 
013 
014 
01s 
016 
017 
018 
019 
020 
021 
022 
023 


025 
026 
027 
028 
029 
030 


-gT- 


Col. 7 
1 


101 DEL=0.0 
Gd TA 103 
102 DEL=.5*XSTEP 
C ....INITIALIZE AS NECESSARY 
103 IST#P=0 
DO 4 I=1,ISIZE 
DO 4 J=1,ISIZE 
4 IGRID(I,J)=0 
NgUT=0 
KOUNT=0 
Gh TA 39 
++++BEGIN PLOTTING SETS MF NI POINTS IN THE MEM@RY GRAPH, IGRID. 
++++IN THIS FIRST STAGE, THE GRAPH IS IN TERMS OF FREQUENCIES. 


a00 


50 GO T(1001,1001,1003) , ITYPE 
1001 Dg 1002 I=1,NI 
1002 ¥(I)=YMIN 
C ...,SET ISTOP Tf 1 IF STOP NUMBER HAS BEEN READ. 
1003 D@ 500 I=1,NI 
IF (X(I)-ST@B) 500,499,500 
499 IS=I-1 
ISTOP=1 
IF(IS)501,62,501 
500 IS=NI 
501 NI=IS 
C ....REJECT POINT IF IT IS QUTSIDE LIMITS OF GRAPH. 
DO 6 I=1,NI 
IF(IDROP-1)71,72,71 
71 IF(((YMIN+(ISIZE-1) *#YSTEP)+. 54YSTEP)-Y(I))7,8,8 
72 IF(((YMIN+(ISIZE-1) #YSTEP)+. 5#YSTEP)-¥(I))7,7,8 


++++FOR FREQUENCY DISTRIBUTIONS, TABULATE IN 1ST INTERVAL ON @RDINATE. 


O31 
032 
033 
034 
035 
036 
037 
038 
039 
040 
041 
042 
043 
044 
045 
046 
047 
048 
049 
050 
051 
052 
053 
054 
055 
056 
057 
058 
059 
060 


-6T- 


Col. 7 
‘ 


7 NOUT=NPUT+L 
Gh Th 6 
8 IF(IDROP-1)91,92,91 
91 IF((YMIN-.5#YSTEP)-Y(I))10,9,9 
92 IF((YMIN-.5*YSTEP)-Y(I))10,10,9 
9 NOUT=NOUT+1 
CH TO 6 
10 IF(ITYPE-2) 210,220,210 
210 IF(IDROP-1) 211,213,211 
211 IF ((XMIN-DEL)-X(I))212,13,13 
212 LF(((XMIN+(ISIZE-1)*XSTEP)+DEL)-X(I))13,14,14 
213 IF ((XMIN-DEL)-X(I)) 214, 214,13 
214 IF(((XMIN+(ISIZE-1) #XSTEP)4+DEL)-X(I))13,13,14 
220 IF(IDR@P-1) 221,223,221 
221 IF((XMIN-XSTEP)-X(I))222,13,13 
222 IF ((XMIN+(ISIZE-1) *XSTEP)-X(I))13,14,14 
223 LE((XMIN-XSTEP)-X(I)) 224,224,13 
224 LF ((XMIN+(ISIZE-1) *XSTEP)-X(I))13,13,14 
13 NOUT=NOUT+L 
Gh TO 6 
C ....L@CATE NEAREST HORIZONTAL COORDINATE. 
14 LOCKX=1 
LoCcY=1 
Dé 16 J=1,ISIZE 
IF (X(1)-( (XMIN+(J-L) AXSTEP)+DEL) ) 16,17 ,9917 
17 IF (IDR@P-1) 16,9917, 16 
9917 LOCX=LOCK+1 
16 CONTINUE 
C ....INCREMENT K@UNT. 
KQUNT=KOUNT+L 


061 
062 
063 
064 
065 
066 
067 
068 
069 
070 
071 
072 
073 
074 
075 
076 
077 
078 
079 
080 
081 
082 
083 
084 
085 
086 
087 
088 
089 
090 


0 


Col. 7 
: 


C ....+LO@CATE NEAREST VERTICAL COPRDINATE. 
Dé 18 J=1,ISIZE 
IF(¥(I)-( (YMIN+(J-1) *YSTEP)+. 5*YSTEP) )18,19,9919 
19 IF(IDR@P-1)18,9919,18 
9919 LOCY=LOcY+1 
18 CONTINUE 
C ....INCREMENT FREQUENCY @F NEAREST GRID INTERSECTION. 
IF (IGRID(LOCY ,LOCX) -32767) 1004, 341,341 
1004 IGRID(LOCY,LOCX)=IGRID (LOCY ,L@CX)+1 
6 CONTINUE 
C ....IF ALL POINTS NOT PROCESSED, RETURN TO MAIN PROGRAM FOR NEXT SET. 
IF(ISTOP)62,5011,62 
5011 CONTINUE 
GO Td 39 
C ....BEGIN REPLACING FREQUENCIES BY PLOTTING SYMBOLS 
C ....FREQUENCY DISTRIBUTIONS AND HIST@GRAMS REQUIRE A SPECIAL PROCEDURE. 
62 GP TM (1009,1009,1300) , ITYPE 
C ....THIS PROCEDURE IS FOR FREQUENCY DISTRIBUTIONS AND HIST@GRAMS 
1009 SAV=0.0 
C ...-LOCATE ZERO POINT ON THE VERTICAL SCALE, IF INCLUDED. 
121 
DP 9991 I=1,ISIZE 
IF(IDROP-1)1014,1015,1014 
1014 IF(YMIN+. 5*YSTEP)9990,9991,9991 
1015 IF(YMIN+.5*YSTEP)9990,9990,9991 
9990 IZ=1Z+1 
9991 CONTINUE 
TQ=0 
C ...,DETERMINE HIGHEST FREQUENCY INCLUDED IN THE VERTICAL SCALE. 
oi 


091 
092 
093 
094 
095 
096 
097 
098 


118 
119 
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Col. 7 
J 


IF ( (YMIN+(ISIZE-1) #YSTEP)+. 5*YSTEP) 39, 39,2101 
2101 FI=I 
IF (IDROP-1) 1016 ,1017,1016 
1016 IF(FI-((YMIN+(ISIZE-1)*YSTEP)+.5¥*YSTEP) )2102,2102,2103 
1017 IF(FI-((YMIN+(ISIZE-1) *YSTEP)+.5#*YSTEP) ) 2102, 2103, 2103 
2102 IT@P=1 
I=I+l 
GO TO 2101 
C ....BEGIN MAIN LOOP FOR FREQUENCY DISTRIBUTIONS AND HISTOGRAMS. 
2103 SAVE=0.0 
DO 1115 J=1,ISIZE 
C ....ACCUMULATE FREQUENCIES OR NOT, ACCORDING TO THE PURPQSE. 
IF (ITYPE-2)1100,1010,1100 
1010 SAV=IGRID(1,J) 
SAVE=SAVE+SAV 
GO TO 1101 
1100 SAVE=IGRID(1,J) 
C ....IS FREQUENCY ABOVE RANGE OF VERTICAL SCALE. 
1101 IF(IDR#P-1)1901,1902,1901 
1901 IF(((YMIN+(ISIZE-1) *YSTEP)+.5*YSTEP)-SAVE)1102,1103,1103 
1902 IF(((YMIN+(ISIZE-1) *YSTEP)+, 5*YSTEP)-SAVE) 1102,1102,1103 
1102 G% TH(1202,1212) ,ITYPE 
1202 IDIF=IGRID(1,J)-ITOP 
KOUNT=KQUNT-IDIF 
NOUT=NOUT+IDIF 
GO TH 1203 
1212 IF(IQ)1203,1213,1203 
1213 NOUT=KOUNT-ITOP 
KOUNT=ITOP 
IQ@1 


121 
122 
‘ZS, 
124 
125 
126 
127 
128 
129 
130 
131 
132 
133 
134 
135 
136 
137 
138 
139 
140 
141 
142 
143 
144 
145 


147 
148 
149 
150 


te- 


Col. 7 
’ 


1203 K=ISIZE-1 
IGRID (ISIZE,J)=ICHAR(1) 
TO TH 1111 

C ....IS FREQUENCY BEL@W THE RANGE OF THE VERTICAL SCALE. 

1103 IF(IDROP-1) 1903 ,1904 ,1903 

1903 IF((YMIN-.5*YSTEP)-SAVE)1106,1104,1104 

1904 IF((YMIN-.5*YSTEP)-SAVE) 1106 ,1106,1104 

1104 GP T6(1220,1240) , ITYPE 

1220 NOUT=NGUT+IGRID(1,J) 
KQUNT=K@UNT-IGRID(1,J) 
GO TH 1250 

1240 IF(J-ISIZE)1250,1241,1250 

1241 NOUT=NOUT+XOUNT 
KQUNT=0 

1250 Dé 1105 I=1,ISIZE 

1105 IGRID(I,J)=ICHAR(100) 
GO TO 1115 

C ....DETERMINE CO@RDINATE @N VERTICAL SCALE. 

1106 K=L 
Dé 1108 I=1,ISIZE 
IF (IDR@P-1) 1905 ,1906,1905 

1905 IF(SAVE-( (YMIN+(I-1) #YSTEP)+, 5*YSTEP))1108,1108,1107 

1906 IF(SAVE-( (YMIN+(I-1) *YSTEP)+, 5*YSTEP)) 1108 ,1107,1107 

1107 K=K+L 

1108 CONTINUE 

C ....PLACE APPROPRIATE SYMBOLS IN COLUMN. 
IF (SAVE) 39 ,9993,1109 

9993 IGRID(K,J)=ICHAR(100) 
IF (ITYPE-2) 9994 ,1109,9994 

1109 IGRID(K,J)=ICHAR(1) 


151 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
178 
179 
180 


Col. 7 


9994 K=K+1 
IF (K-ISIZE) 9992,9992,1908 
9992 D@ 1110 I=K,ISIZE 
1110 IGRID(I,J)=ICHAR(100) 
1908 K=K-2 
1111 IF(K-1Z)1115,1909,1909 
1909 Dé 1112 I=1Z,K 
1112 IGRID(I,J)=ICHAR(2) 
IF(IZ-1)1115,1115,1911 
1911 K=1Z-1 
Dé 1912 I=1,K 
1912 IGRID(I,J)=ICHAR(100) 
1115 CONTINUE 
GH TA 1500 
C ....THIS PROCEDURE IS FOR ITYPE NUMBER 3. 
1300 1Q=0 
DO 20 I=1,ISIZE 
Dé 20 J=1,ISIZE 
K=IGRID(I,J) 
IF(K)22,21,22 
21 IGRID(1,J)=ICHAR(100) 
GH TH 20 
22 IF(K-100)24,23,23 
23 IGRID(I,J)=ICHAR(101) 
C ....PRINT THE EXACT FREQUENCY WHEN @VER 99. 
IF(IQ)100,34,100 
34 WRITE (19,36) ICHAR(101) 
36 FORMAT(1H1,59HIN THE GRAPH BEL@W, FREQUENCIES QVER 99 ARE REPRESEN 
ITED BY ,A2,1H.,/,38H A LIST OF THESE FREQUENCIES FOLLOWS. ,//,25X, 
218X, 20X, 1HY,15X,9HFREQUENCY //) 


181 
182 
183 
184 
185 
186 
187 
188 
189 
190 
191 
192 
193 
194 
195 
196 
197 
198 
199 
200 
201 
202 
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Col. 7 
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IQ=1 
100 XSCAL=XMIN+(J-1) *XSTEP 
YSCAL=YMIN+(J-1) *YSTEP 
WRITE(1I@,37)XSCAL,YSCAL,K 
37 FORMAT(LH ,20X,2(F10.4,10X),15) 
Gd TH 20 
24 IGRID(I,J)=ICHAR(K) 
20 CONTINUE 
C ....BEGIN PRINTING GRAPH(S). 
1500 Dé 390 NPRNT=1,NCOPY 
C ....PRINT GRAPH. 
WRITE(IO,25)YTITL 
25 FORMAT(1H1,//,1X,18A4) 
L=4 
DO 26 K=1,ISIZE 
I=(ISIZE+1)-K 
IF(L) 28, 27,28 
27 YSCAL=YMIN+(I-1) *YSTEP 
WRITE (1,29) YSCAL, (IGRID(I,J) ,J=1,ISIZE) 
L=5 
Gh TO 26 
28 WRITE(16,30) (IGRID(I,J) ,J=1,ISIZE) 
26 L=L-1 
29 FORMAT(1H ,7X,F10.4,2H +,50A2) 
30 FORMAT(1H ,17X,2H +,50A2) 
IF (ISIZE-50) 605,604,605 
604 WRITE(I9,31) 
31 FORMAT(1H ,18X,101H+ +- ---+----+----+----+-- 
l--+----+----+----+---- 
Gd Td 606 





211 
212 
213 
214 
215 
216 
217 
218 
219 
220 
221 
222 
223 
224 
225 
226 
227 
228 


230 
231 
232 
233 
234 
235 
236 
237 
238 
239 
240 
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605 WRITE(I9,311) 
311 FORMAT(1H ,18X, 81H+ +----+----+----+----+-- 
1--+----+----+----) 
606 IQ=(ISIZE/5)-1 
Dé 5 I=1,1Q 
5 X(1)=XMIN+I*5*XSTEP F 
WRITE (19,32) KOUNT ,XMIN, (X(I) ,I=1,1Q) 
32 FORMAT(1H ,15,10X,10F10.4) 
WRITE(I9,33)XTITL 
33 FORMAT(1H ,7HENTRIES,//,21X,18A4) 
IF(NOUT)390,390,38 
++++ERR@R MESSAGE FOR THE UNLIKELY EVENT @F A FREQUENCY OVER 32766. 
341 WRITE(I6,3401) 
3401 FORMAT(1HO,65HFREQUENCY OF AN INTERSECTION EXCEEDS 32766, EXECUTIO 
IN TERMINATED. ) 
Gh TO 39 
«+++PRINT NUMBER OF POINTS REJECTED DUE T@ BEING BEYOND SCALE LIMITS. 
38 WRITE(I6,40)NOUT 
40 FORMAT(1H1,//,15,34H POINT(S) BEYOND THE SCALE LIMITS,/,51H AND TH 
1EREFQRE TRUNCATED QR OMITTED FROM THE GRAPH.) 
390 CONTINUE 
39 IC(1)=1Gp 
IC(2)=NI 
IC(3)=ITYPE 
Ic(4)=10 
IC(5)=ISIZE 
1c (6)=NCOPY 
IC (7)=ISTOP 
IC (8)=IDROP 
IC (9)=NgUT 
IC(10)=KdUNT 
RETURN 
END 


241 
242 
243 
244 
245 
246 
247 
248 
249 
250 
251 
252 
253 
254 
255 
256 
257 
258 
259 
260 
261 
262 
263 
264 
265 
266 
267 
268 
269 
270 
271 
272 
273 
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VIIL. EXAMPLE #1: X-Y PLOTTING 


Following is the first of five complete examples which 
should aid the user in applying the program to his data. Each 
example illustrates a different type of application. 


Assumed data 


Suppose that we have 22 observed pairs of values, the first 
member of a pair being the value of a variable X and the second 
member the value of a variable Y. Each observed pair of values is 
punched on a data card, with the X-value in columns 5 through 10 
and the Y-value in columns 26 through 30. Although no decimal 
points are actually punched on the cards, columns 10 and 30 of each 
card are implied decimal places. Thus, the digits 25 punched in 
columns 9 and 10 of a card represent an X-value of 2.5 . The 
observed values and their positioning on the cards are shown in 
Figure 5. 


Purpose 


Our purpose in this example is to obtain a simple X-Y plot 
of the data, with the variable X on the horizontal scale and the 
variable Y on the vertical scale. The symbol plotted at an inter- 
section of grid lines should be the number of points for which this 
is the nearest intersection. If the number of points falling at an 
intersection exceeds 99 (which is an impossibility in this case 
since we only have 22 points), an asterisk should be placed at that 
intersection. 


We choose zero as the midpoint of the lowest scale interval 
on both axes. We choose .5 as the scale interval on the horizontal 
axis and 1.0 as the scale interval on the vertical axis. We wish the 
resulting graph to fit on an 83 x ll-inch page, so we choose the 
smaller size. Our example is now constructed so that the first data 
point and the last three data points fall outside the scale limits and, 
therefore, will be omitted from the graph. The code number of the 
printer is 3 since this example is to be run on an IBM 1130 system 
with an IBM 1132 printer. Finally, we do not require any additional 
copies of the graph. 


Main program 


The main program appropriate for this example is shown in 
Figure 6. With respect to the READ card, note that the code number 
of the IBM 1442 card read punch is 2 . 


Control cards 
The completed control cards for this example are shown in 


Figure 7. We use 9999.9 as the flag number since this number does 
not occur among our observed values for the variable X. 


my 
2 
d SSSRSASLARNS RA AASSRSISSSANSRRSSSSSLRLLS 
250 500 
0 0 
bi 10 
10 20 
10 20 
20 40 
40 80 
50 100 
50 100 
50 100 
70 140 
80 160 
100 200 
100 200 
110 220 
120 240 
130 260 
135 270 
145 290 
200 400 
210 420 
200 900 
Figure 5, Assumed data for example #1. Arrows at bottom and top of figure show the 


location of the implied decimal points. 
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Results 


The results of running example #1 are shown in Figures 8a 
and 8b. The first page contains the graph, and the second page 
contains supporting information regarding the number of points 
which were omitted because they were beyond the scale limits. 


Column 


number 


Figure 6. 


Aron 
oooor 


12 
14 
16 
18 
20 
22 
24 
26 
28 
30 


DIMENSION A(18) ,B(18) ,X(20) ,¥(20) ,IA(50,50) ,1B(101) , 1¢(10) 
Te(1)=1 

GO TA 4 

IF(1C(7)-1)2,5,2 

NI=IC (2) 

READ (2,3) (X(I) ,¥(I) ,I=1,NI) 

FORMAT (5X,F5.1,15X,F5.1) 

CALL PLOT(A,B,X,Y,C,D,£,F,G,H, IA, 1B, IC) 
Ic(1)=2 

Gp TA 1 

CALL EXIT 

END 
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Main program for example #1. 


Column 
number 


Control 
card 


SYMBOLA 
SYMBOLB 
SYMBOLC 
XTITLE 
YTITLE 
OPTION 


STOP 


Figure 7. 


123.4 5 6 7 8 910111213141516171819202122232425 2627 28293031323334353637 383940 
4142434445464748495051525354555657585960616 26364656667 68697071727374757677787980 


81828384858687888990919293949596979899 * 


X-VALUE 
Y-VALUE 
0.0 0.5 0.0 1.0 9999.9 2 ee ee es | at 
99999 


Completed control cards for example #1. 
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Y=VALUE 


35.0000 


3020000 


“eg aan3ty 


2520000 


2020000 


1540000 


“T# eTduexe roxy yderd yndqng 
$9447.44 $F £44444 FHF H OHHH 44H HH 4H H+ HHH 44+ + HF 


10.0000 3 
1 
520000 
1 
2 
1 

020000 1 

| iat Siete Sie, Si So A a 

18 - 0.0000 245000 520000 725000 10.0000 1225000 15.0000 1725000 
ENTRIES 


X=VALUE 
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& POINT(S) BEYOND THE SCALE LIMITSs 
AND THEREFORE TRUNCATED OR OMITTED FROM THE GRAPHe 


Figure 8b. Supporting information printed for example #1. 
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IX. EXAMPLE #2: HISTOGRAM 
Assumed data 


Suppose that we have 75 observations on a single variable. 
The measurements are only to the nearest unit, so that there are no 
decimal places. Five observed values are punched on each card; each 
observed value being allotted a field of five columns in which it is 
right justified. The first value on a card occupies colums | 
through 5, and five blank columns separate adjacent fields. The 
data and the layout on the cards are shown in Figure 9. 


Purpose 


Our purpose is to construct a’ histogram of these data. A 
set of two dashes (minus signs) should mark the top of a column, and 
the column should be shaded below with sets of two dots (decimal 
points or periods). We choose 0. as the midpoint of the lowest 
interval on both scales, and 1.0 as the scale interval on both scales. 
Again, we want the smaller size of graph, and do not require any 
duplicate printouts. As before, the code number of the line printer 
on the IBM 1130 systemis 3 . 


Main program 


The main program appropriate for this case is shown in 
Figure 10, It is only necessary for the program to read X-values 
since we are dealing with a single variable which is to appear on the 
horizontal axis; however, the "Y(20)" must still appear on the 
DIMENSION card. 


Control cards 


The coding for the control cards is illustrated in Figure ll. 
Blank cards are simply inserted in the places of the SYMBOLB and 
SYMBOLC control cards. 


Results 


The histogram produced by the computer for this example 
run is shown in Figure 12. In this case, the zero line on the 
vertical axis is superfluous and might better be deleted by using 
1.0 as the midpoint of the first interval on the frequency scale. 
With frequency scale intervals greater than one, it is advisable 
to include zero in the first interval. A zero frequency (or no 
observations for a given interval) will then be indicated by a 
blank at the corresponding position in this first line. 


Column 
number 


SSSSRAASASAARARAASASS 
25 25 25 25 
25 25 25 25 
25 25 25 25 
25 25 25 25 
10 10 10 10 
15 15 15 15 

2 if 2 2 

3 4 4 5 
20 20 20 20 
20 20 20 20 
30 30 30 30 
30 30 30 35 
35 35 35 35 
35 35 35 35 
35 35 35 35 


Figure 9, Assumed data for examples #2 and #3. 
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Column 


number 


NK 


Figure 10. 


DIMENSI@N A(18) ,B(18) ,X(20) ,¥(20) ,TA(50,50) ,IB(101) ,10(10) 
IC(1)=1 

GO Td 4 

IF(IC(7)-1)2,5,2 

NI=1C (2) 

READ (2,3) (X(1) ,I=1,NI) 

FORMAT (5(F5.0,5X)) 

CALL PLOT(A,B,X,Y,C,D,E,F,G,H,IA, 1B, 1C) 
Ic(1)=2 

GO TA 1 

CALL EXIT 

END 
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Main program for examples #2 and #3. 


azequnu 
uunto9 


Control 


card 


SYMBOLA 


SYMBOLB 


SYMBOLC 


@BSERVED VALUE 


XTITLE 


a 


FREQUENCY 


YTITLE 


1.0 0.0 1.0 99999. os oS Ss 1,4 


0.0 


OPTION 


99999 


STOP 


Control cards for example #2. 


Figure 11. 


“ZT eandyy 


*Z# eTdwexe 103 aaqndwos ay fq paonpoad wea8oasty 


75 
ENTRIES 


FREQUENCY 


3520000 


3020000 


25420000 


20.0000 


1520000 


1020000 


520000 


020000 


HHH EHP HEHEHE HEHEHE HS Hee eee tHe eee eeetee tee teteeee tee + 


re 





oe 
ee 
ee 
oe 
oe 
oe 
oe 
oe 
oe 
es 
oe 
oe 
oe 
oe 
oe 
oe 
ee 
oe 
oe 
oe 


oe 
ee 
oe 
ee 
ee 
ee 
ee 
oe 
oe 
ee 
oe 
oe 
ee 
oe 
oe 
oe 
ee 
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ee ee ee ee 


0.0000 


5.0000 


10.0000 


15.0000 


20.0090 


OBSERVED VALUE 


25.0000 


30.0909 


3529900 
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X. EXAMPLE #3: CUMULATIVE FREQUENCY DISTRIBUTION 
Assumed data 


We use the same data as for example #2 (Figure 9). 


Purpose 


Our purpose is to construct a cumulative frequency distribu- 
tion by using the data from example #2. We will use an asterisk (*) 
as our plotting symbol, and there will be no shading below. Also, 
we shall use 2.0 as the scale interval on the vertical axis 
(frequency scale). All other specifications remain the same as for 
example #2. 


Main program 


No change is required in the main program (Figure 10). 
Control cards 


The coding for the control cards is shown in Figure 13. 
Again, the SYMBOLB and SYMBOLC control cards are blank. 


Results 


The results of running example #3 on an IBM 1130 computer 
are shown in Figure 14. We note that using a frequency scale inter- 
val of 2.0 and leaving column 53 of the OPTION control card blank 
effectively makes the numbers printed on this scale lower limits of 
the frequency classes. This is the result of including the upper 
limit of the class and excluding the lower limit when the numbers can 
assume only integral values. 


Column 
number 


Control 
card 


SYMBOLA 
SYMBOLB 
SYMBOLC 
XTITLE 
YTITLE 
OPTION 


STOP 


Figure 13. 


OBSERVED VALUE 


-0%- 


CUMULATIVE FREQUENCY 
0.0 1.0 0.0 2.0 99999, me 222 2 


99999 


Coding of control cards for example #3. 


"€# eTdwexe a0zy azagnduos ayy Aq 
paonpoad uozjnqzaqstp Aousnbeay aataeqnwng ‘yy eanyy 


a> 
ENTRIES 


CUMULATIVE FREQUENCY 


70.0000 


6020000 


5020000 


4020000 


30.0000 


2020000 


1020000 


0.0000 


+ 


tHe ee Hotere eeeeeeeeeeeeeeeeeee ee eer ee ete 


+ 


a rr 


0.0000 


5.20000 


* 2 * 


se eee 


ee ee 


eee HH 


ee RH He 


10.0000 


15.0000 20,0000 


OBSERVED VALUE 


25.0000 


30.0000 


eee RK 


ae + 
rs 
1 
35.0000 
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XI. EXAMPLE #4: CONTINGENCY TABLE 
Assumed data 


Assume that 20 subjects have been classified with respect 
to two attributes, which we call attribute A and attribute B. Assume 
further that there are two categories for each attribute, coded 
respectively as 1 and 2 . The classification data for each subject 
have been punched on a card with the code for attribute B in column 5 
and the code for attribute A in column 10.- These hypothetical data 
cards are illustrated in Figure 15. 


Purpose 


Our purpose is to fill in the cells of a 2 x 2 contingency 
table corresponding to these data. Attribute A is to be arrayed along 
the left margin and, therefore, corresponds to an ordinate or Y-value, 
Attribute B is to be arrayed along the lower margin and, thus, cor- 
responds to an abscissa or X-value. The number 1 is appropriate for 
the midpoint of the lowest scale interval along both margins, since 
this is the code for the first category of both attributes. We want 
the code number 2 to be printed below (or beside) the frequency for 
the second category of both attributes. Since the program prints 
only every fifth scale value, we should be able to achieve this by 
using a scale interval of .2 for both margins. The program will give 
the cell frequencies and the total frequency, but the other marginal 
frequencies must be calculated manually from the cell frequencies. 
Also, frequencies over 99 cannot be printed in the table, but a list 
of such frequencies (if any) would precede the first copy of the 
table. 


Main program 


The main program appropriate for example #4 is shown in 
Figure 16. 


Control cards 

The coding for the control cards is illustrated in Figure 17. 
Note that these control cards would cause any zero frequencies in 
the table to appear as blanks. 


Results 


The output of an IBM 1130 computer, for example run #4, is 
shown in Figure 18. 


“CT ean3tq 


*yg eTdwexs 10y elep pounssy 


FPREN NNN NEF NNNNNNNE RENE 


NPR PR NER RP RRP RR NER ERNE 
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Column 
number 


Column 
number 


ssesoruseccaseaonsea 
SCOSCCHH AA AANANNANMAMMOM 


40 
42 
44 
46 
48 
50 


DIMENSION A(18) ,B(18) ,X(20) ,Y(20) ,1A(50,50) ,1B(101) ,1C(10) 
Ic(1)=1 
GO TO 4 
1 IF(IC(7)-1)2,5,2 
2 NI=ICc(2) 
READ (2,3) (X(I) ,¥(I) ,I=1,NI) 
3 FORMAT (4X,F1.0,4X,F1.0) 
4 CALL PLOT(A,B,X,Y,C,D,E,F,G,H,IA,1IB,1C) 
Ic(1)=2 
GO TH 1 
5 CALL EXIT 
END 


py 


Figure 16, Main program for example #4. 


Column 
number 


Control 
card 


SYMBOLA 
SYMBOLB 
SYMBOLC 
XTITLE 
YTITLE 
OPTION 


STOP 


Figure 17, 


123.45 6 7 8 910111213141516171819202122232425262728293031323334353637383940 
4142434445464748495051525354555657 585960616 2636465666 768697071727374757677787980 


81828384858687888990919293949596979899 ., 


ATTRIBUTE B 
ATTRIBUTE A 
1.0 0.2 1.0 0.2 9. Ls) Oe a A 


Coding of control cards for example #4. 
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ATTRIBUTE A 


840000 


720000 


“ST ean3Ta 


600000 


500000 


420000 


320000 


“yy eTdwexe tof 1eqyndmos O¢TT WAI JO 3ndqno 


220000 


1.0000 


20 
ENTRIES 


CEASE StS € HSS SSF OHSS EH AT SS FH tS 44 FH tS t+ HS 
-94- 


2 5 
5 8 

c+ = = Se SS Sy eS SS he = ee Se Ry See SS fee ee hp ee 
1.0000 220000 320000 40000 520000 60000 720000 8.0000 


ATTRIBUTE B 
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XIL. EXAMPLE #5: PLOTTING A FUNCTION 
Assumed data 


No input data are required. 


Purpose 


Our purpose is to plot the area of a circle in square feet 
as a function of the diameter in inches, The midpoint of the first 
interval on both scales will be 0.0 . The scale interval for dia- 
meter will be .5 inch. The scale interval for area will be .1 
square foot. An asterisk (*) is the symbol to be used in plotting 
the function. 


Main program 


The appropriate main program is shown in Figure 19, 
Control cards 


The coding for the control cards is shown in Figure 20. 
The SYMBOLB and SYMBOLC control cards are blank. In this case, the 
STOP control card is omitted entirely since the number in columns 
41 through 50 of the OPTION card is sufficient to signal the end of 
the operation, 


Results 
The graph produced by the computer for example #5 is shown 


in Figure 21. The apparent discontinuities are due to plotting a 
continuous function in discrete intervals. 


number 


Column 


Sseecissegaseaesnesesoseessszesssiseseesees 
SSSCOAHAAANNNNNMMAMNMNANSSTSTIANANNOSOOCORARARARD 


DIMENSION A(18) ,B(18) ,X(20) ,¥(20) ,IA(50,50) ,1B(101) ,1C(10) 
Ic(1)=L 
N=0 
1 CALL PLOT(A,B,X,Y,C,D,E,F,G,H,IA,IB,IC) 
Ic(1)=2 
IF(IC(7)-1)2,6,2 
2 Dé 3 I=1,10 
X(1)=C+NaD 
¥(I)=.005454154#X (I) *X (I) 
3 N=Nt+1 
IF(N-IC(5))5,5,4 
X(1)=¢ 
GO TH 1 
CALL EXIT 
END 


aut 


Figure 19. Main program for example #5. 
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Column 
number 


Control 
card 


SYMBOLA 
SYMBOLB 
SYMBOLC 
XTITLE 
YTITLE 


OPTION 


Figure 20. 


SS SBS SARA RS RR RAS ASLISSRARRASSSSSLRELLS 
* 
DIAMETER @F CIRCLE (INCHES) 
AREA (SQ. FT.) 5 
1 
0.0 0.5 0.0 0.1 0. Jd) 2d a 


Coding of control cards for example #5. Note that the STOP control card is omitted, 


AREA (SQe 


345000 

& 340000 

Cc 

n 

e 

r 
25000 

§ 

a 

3 

ec 

i 

8 

g 

9 220000 

5 

i=] 

e 

is 

rn 

8 

a 145000 

ad 

& 

a 

sc 

2 

® 

a 140000 
025000 
0.0000 

40 

ENTRIES 


FTe) 


SHEE HE HEHEHE HEHEHE HEHE HHH eee eee eee eee eeereeette 
* 
=ge= 


eee 
+e ee eS 


a 


0.20000 25000 5.0000 7.5000 10,000¢ 12.5000 520000 1725000 


MTAMFTFR OF CIRCLE Lt INCHFS) 
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XTIL. STORAGE REQUIREMENTS, TIMING, AND PRECISION 


All information given in this section applies to the [BM 
1130 version IT monitor system using the one word integers and 
extended precision features of the FORTRAN compiler. 


Core storage requirements (in terms of 16-bit words) for 
SUBROUTINE PLOT are as follows: variables = 52, program = 2372. 
All of- the example main programs require 2862 words for variables, 
and 130 words or less for program. Storage requirements for 
another computing system may differ considerably from these figures 
because the one word integers feature of the IBM 1130 FORTRAN 
compiler reduces the storage requirements, whereas the necessity of 
storing the set of subroutines for real arithmetic increases the 
requirements. 


The IBM 1130 system did not provide a printed record of 
processing times. Furthermore, with an IBM 1132 printer, compila- 
tion time is highly dependent on whether or not a listing is 
required. After SUBR(UTINE PLOT had been stored in object form on 
disk, the time required to compile a main program and process an 
example ranged from 2 to 3 minutes. Approximately | minute of this 
was taken to load the object form of the subroutine into core from 
the disk. These times may seem long in comparison to IBM 360 
system speeds, but the cost of the IBM 1130 is proportionately low. 


Especially with an IBM 1130 computing system, it is advis- 
able always to use extended precision for real variables. On such 
a system, extended precision provides three 16-bit words for each 
real variable, thus increasing the number of significant digits 
from 7 to 10. 
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A FORMAT is 
describes the type, 
punched, or printed 
card is as follows: 

Columns 


4 


2 through 5 


7 through 13 


14 through 72 


73 through 80 
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APPENDIX: SIMPLE FORTRAN PORMATS 


a non-executable FORTRAN program statement which 
form, and positioning of data to be read, 
by a computer. The general layout of a FORMAT 


Information 


Leave this first column of the FORMAT card 
blank. For the sake of completeness, it should 
be noted that the FORMAT identification number 
described next is allowed to extend into column 
i under most circumstances, To avoid considera~ 
tion of the special conditions which can arise, 
however, it is preferable to avoid using this 
column. 


This is an integer number, i.e,, no decimal 
point, which identifies the particular F@RMAT 
statement. Leading zeros do not have to be 
punched, and the positioning of the number is 
not critical. Thus, a two-digit number could 
eccupy columns 2 and 3, 3 and 4, or 4 and 5. 
However, the neatness and readability will be 
enhanced if all FORMAT statements in a program 
have the last digit of their identifying number 
in the same column. 


The number in this column is called a continua- 
tion code number. Column 6 of the first card in 
a FORMAT statement may either contain a zero or 
be left blank. Column 6 of any additional cards 
must contain one of the digits 1 through 9. The 
order of numbers in column 6 of cards within a 
FORMAT statement does not necessarily have to 
correspond with the position of the cards in the 
statement; however, interpretation is simplified 
if this convention is adopted. 


FORMAT ( 


These columns contain field descriptions, commas, 
slashes, or parentheses. These are the working 
elements of the FORMAT statement. Their uses 
are explained below. The FORMAT statement must 
always end with a right parenthesis in one of 
these columns. 


The computer completely disregards any informa— 
tion which is placed in these columns. 
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The term "field" designates a set of adjacent columns on 


the input or output medium, whether it be ecards or the printed page 
A generalized field description can be represented as follows: 


where: 


rne.d 


x (meaning "repeat") is the number of times the current 
field description is to be used before the next one is 
considered. 


n is a code letter designating the type of information 
which the field contains. This code letter is the primary 
component of the field description, the other components 
serving as modifiers. Certain ones of the modifiers are 
omitted, depending on the code letter; however, a code 
ietter must always be present in a field description. 


¢ (meaning "columns") is the total number of columns in 
the field, including an allowance for sign and/or decimal 
point if required. 


d (meaning "decimal places") is the number of decimal 
places in the field, not counting the decimal point itself. 
Whenever the d modifier is present in a field description, 
it must be preceded by a decimal point as indicated in the 
general form above. 


The field descriptions are separated from each other by commas. 


Several 
however 
five: 


code letters for different types of fields are available; 


» one can usually operate satisfactorily with the following 
I, F, A, H, and X. The details regarding purposes and uses 


of these codes are as follows: 


Code Purpose 
I This code signifies that the field contains an 


integer number. An integer number cannot con- 
tain any fractional part or decimal point, there- 
fore, the d modifier and its preceding decimal 
point are always omitted from an I-field descrip- 
tion. In addition, the r modifier may be omitted 
if the field description is not to be repeated 
before the next description is considered. 


This code indicates that the field contains a 
number which either does or could contain a 
decimal point. The F-field description requires 
all modifiers, except that the r modifier may 
be omitted if the field description is not to be 


Code 





F 


(cont 'n) 


Purpose 


repeated before the next one is considered. If 
data are being printed or punched, one column 
will be taken up by a decimal point with d 
decimal places following the decimal point. The 
decimal places will always appear as far to the 
right as possible in the field. A number so 
positioned is said to be right justified. When 
data are being read from cards, the d modifier 
covers the case where no decimal point is 
punched but the last d columns in the field are 
to be interpreted as decimal places. A decimal 
point actually punched anywhere in the field on 
the data card will completely override the d 
modifier. 


This code signifies that the field contains 
alphameric information, Letters, numbers, and 
punctuation marks are all included under the term 
alphameric; however, numbers read under this 
designation cannot be used in arithmetic opera- 
tions, This code is primarily useful for 
manipulating names and titles which are to be read 
from data cards. The d modifier and its preceding 
decimal point are always omitted from A-field 
descriptions. Unless one is familiar with the 
particular computer being used, it is safest to 
make the c modifier no larger than 2 . This may 
frequently be wasteful of computer storage, but 
should avoid the possible omisSion of some 
characters. The r modifier may be omitted if the 
field description is not to be repeated. 


The H code is used to specify carriage control 
(explained later) and to transmit alphameric 
information to or from the FORMAT statement 
itself. The H code is basically a single column 
operator, and, therefore, both c and d modifiers 
with their intervening decimal point are always 
omitted. On input, the characters in the next 
x columns will be inserted into the FORMAT 
statement immediately after the H. On output, 
except for carriage control, the r characters 
immediately following the H will be punched or 
printed just as they appear in FORMAT statement. 


This code indicates that a column on the input 
or output medium is to be skipped. The X code 
is basically a single column operator and, there- 
fore, both the c and d modifiers with their 


Code Pucpose 


x (cont'n) intervening decimal point are always omitted. 
The r modifier indicates the number of columns 
that are to be skipped before the next field 
description is considered. 


A / is inserted in the FORMAT statement whenever the com- 
puter is to skip the remainder of the current line or card and 
proceed to the first column of a new Line or card. Commas before 
and after a / are optional. If the FORMAT controls a printing 
operation, a carriage control character must be provided for the 
new line. It also is possible to repeat a group of field descrip- 
tions by enclosing the group in parentheses and placing an r 
modifier before the left parenthesis. 


Up to this point, carriage control characters have been 
mentioned without explanation. If output is the printed page, the 
computer must know where the line is to be placed in relation to the 
last line printed. To specify this, an H code with r modifier equal 
to 1 is placed at the beginning of the FORMAT statement and after 
each transfer to a new line. The character following the H is called 
a "carriage control character" and determines the spacing between 
lines. If the number | follows the H, printing will begin at the top 
of a new page. The number 0 causes a double space, and a blank 
causes a single space. Other carriage control characters are avail- 
able, but usually will not be needed. 


In operation, a FMRMAT statement functions as follows. Each 
input or output statement in the program specifies a list of items 
to be transmitted and the number of a FORMAT statement which is to 
control the transmission of the information. All columns of the 
cards or lines involved must be accounted for in the FORMAT statement 
up to the point where transmission stops. The computer starts at the 
left of the next available card or line of the input or output medium. 
The field descriptions in the FORMAT statement are used in order from 
left to right until no further items of information remain to be 
transmitted. If the last field description is used before trans— 
mission is completed, the computer goes to a new line or card on the 
input or output medium and then returns to the nearest left paren- 
thesis in the FMRMAT statement. The field descriptions after this 
parenthesis are again used in order from left to right until all 
items in the list have been transmitted or the last field description 
is again encountered. In the latter case, the above process is again 
repeated, etc. 


A complete understanding of FORMAT statements requires some 
knowledge of FORTRAN input/output statements and variable types, which 
are not discussed here. However, a common problem is to apply an 
existing program to data which meet all of the requirements of the 
program except that they are not positioned properly on the data cards. 
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If the new data occur in the same order on the cards as the original 
data, it is only necessary to alter the appropriate FORMAT state 
ments by changing the modifiers of the I- and F~ field descriptions 
and inserting or removing X codes or / 's as required, If the 
items of data do not occur in the proper order, it also is necessary 
to reorder their names in the READ statements and their correspond- 
ing field descriptions in the FORMAT statements. 


For example, suppose that the relevant statements in the 
original program are as follows: 


(Col. 7) 


READ(3,1) (JF(L) ,2(1) ,Y(L) , I=L, NT) 
1 FORMAT (15,5X,FL0.5,10X,F10.5) 


This set of statements specifies that three items of data are to 
be read from each data card, The first item, named JF, can never 
have any fractional part; for example, it might be a frequency. 
This datum occupies columns | through 5 of each card, Columns 6 
through 10 of the card do not contain relevant information. The 
second datum, named Z, occupies columns 11 through 20 of the card. 
If no decimal point is punched, columns 16 through 20 are to be 
interpreted as decimal places. Columns 21 through 30 do not contain 
relevant information. The third datum, named Y, occupies columns 
31 through 40 and is similar in form to the second. The READ 
statement given above also specifies that NI cards will be read 
each time the statement takes effect. We do not need to worry 
about the actual value of NI since it is determined elsewhere in 
the program. 


Now suppose that we wish to adapt this program to a new set 
of data which has the following characteristics. Hach value of the 
variable represented by JF now occupies columns 8 through LO and is 
the only item on the card. The values of Z and Y which would have 
appeared on the card under the original scheme now appear on a 
second card. The value represented by Z appears in columns 6 through 
10, with columns 9 and 10 being decimal places. Columns 11 through 
20 do not contain relevant information, Finally, the value of Y 
appears in columns 21 through 30. This value of Y is recorded only 
to the nearest unit so that there are no decimal places. Since the 
items of data occur in the same order as before, the READ statement 
does not require revision. The appropriate revision of the FORMAT 
statement is: 


(Col. 7) 
' 


1 FORMAT (7X,13/5X,F5.2,10X,F5.0) 
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Under this new scheme, NI sets of two cards will be read each time 
the READ statement takes effect. The value of JF for a given case 
will be taken from the first card, and the corresponding values of 
Zand Y will be taken from the second card. Note that the F-field 
designation is retained for the variable Y even though no decimal 
places are recorded. The program will not function properly unless 
the original code letters are retained. 


As a further illustration, suppose that the new data set 
differs from that just described only in having the positions of 
the Z and Y fields interchanged, so that Y precedes Z on the 
second card. The corresponding interchanges required in the READ 
and FORMAT statements are as follows: 


(Col. 7) 
' 


READ (3,1) (JF(1) ,Y(1) ,Z(1) , I=1,NI) 


1 FORMAT(7X,13/5X,F5.0,10X,F5.2) 


